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Abstract The Helioseismic and Magnetic Imager (HMI) onboard the Solar 
Dynamics Observatory (SDO) is designed to study oscihations and the mag- 
netic field in the solar photosphere. It observes the full solar disk in the Fe I 
absorption line at 6173 A. We use the output of a high-resolution 3D, time- 
dependent, radiation-hydrodynamic simulation based on the CO^BOLD code to 
calculate profiles F(A,x,?/, t) for the Fe I 6173 A line. The emerging profiles 
F(A,x,7/,t) are multiplied by a representative set of HMI filter transmission 
profiles i?i(A, 1 < i < 6) and filtergrams Ii{x,y^t;l < i < 6) are constructed 
for six wavelengths. Doppler velocities Vhmi(^,^,^) are determined from these 
filtergrams using a simplified version of the HMI pipeline. The Doppler velocities 
are correlated with the original velocities in the simulated atmosphere. The cross- 
correlation peaks near 100 km, suggesting that the HMI Doppler velocity signal 
is formed rather low in the solar atmosphere. The same analysis is performed for 
the SOHO/MDI Ni l line at 6768 A. The MDI Doppler signal is formed slightly 
higher at around 125 km. Taking into account the limited spatial resolution of the 
instruments, the apparent formation height of both the HMI and MDI Doppler 
signal increases by 40 to 50 km. We also study how uncertainties in the HMI 
filter-transmission profiles affect the calculated velocities. 
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1. Introduction 



The Helioseismic and Magnetic Imager (HMI: Schou et al. 2011) onboard the 
Solar Dynamics Observatory observes the full solar disk in the Fe I absorption 
line at 6173 A to measure oscillations and the magnetic field in the solar pho- 
tosphere. The oscillation measurements are used for helioseismic studies of the 
solar interior, the magnetic-field measurements for studies of solar activity. 
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HMI supersedes the Michelson Doppler Imager (MDI: [Scherrer et al 1995 ) 



onboard the Solar and Heliospheric Observatory (S OHO: [Domingo et aL\\[996 )^ 
which had a similar design and observing strategy. HMI offers significantly im- 
proved spatial resolution (4k x 4k detectors vs. a Ik x Ik detector on MDI), 
increased image cadence (45 vs. 60 seconds on MDI), full Stokes capabilities 
(MDI measured only the longitudinal component, or Stokes V), and all observ- 
ables are transmitted to the ground without any further onboard processing, 
which allows for a much better correction of instrumental effects than would be 
possible with more limited on-board resources. Another significant change is the 
use of a different spectral line. While MDI observed the Ni I 6768 A line, which 
is also used by the GONG project, for HMI the Fe I line at 6173 A was selected. 
The reason for this switch is that the Fe I line is magnetically more sensitive and 
offers better performance for vector magnetic-field measurements than the Ni I 



line (Norton et a/., 2006). 


The Fe I 6173 line has 


□een used extensively in solar magnetometry (e.^.'Baur| 
|2000, Bello Gonzalez et al.^ 2008, 2009). Its formation 


et al. 1980 Muller et al. 


process has been studied 


by e.g. Bruls et al. (1990), Norton et al. (2006), and 


Bello Gonzalez et al. (2009). Bruls et al. (1990) quote 276 km and 238 km for the 



formation height of the line core in the VAL-C atmosphere in LTE and NLTE, 
respectively. Norton et al. (2006) found a slightly higher value of 302 km for 
the NLTE formation height of the line core in their calculations, which were 
also performed for the VAL-C atmosphere. Bello Gonzalez et al. (2009) present 
velocity response functions, which show a broad maximum between 200 and 300 
km for the line core and a maximum around 120 km for the center of gravity of 
the Fe 6173 line. 

The HMI Doppler signal is not a measure of the Doppler core velocity. It is 
determined from a combination of six filtergrams taken around the Fe 6173 A 
line. An atlas profile of the Fe 6173 line together with representative transmission 
profiles of the six HMI filters is shown in Figure [T] The nominal HMI filter 
bandwith is 76 mA and the nominal tuning range 680 mA. This is accomplished 
with a set of successively narrower bandpass filters, consisting of a front win- 
dow, a broadband blocking filter, a five-element Lyot filter and two Michelson 
interferometers. 

The aim of this study is to determine the formation height of the HMI Doppler 
signal, which is important for many science questions. Spectral lines are of course 
not formed in a narrow, static layer, but over a broad, highly variable range as 
the overshooting convection and p modes push the atmosphere around. The term 
"formation height" should therefore be understood as a broad average, not as a 
sharply defined geometrical height in a ID, plane-parallel atmosphere. 



2. 3D Model Atmosphere and Line Synthesis 

We study the formation of the HMI Doppler signal with the help of 3D time- 
dependent radiation-hydrodynamic simulations which provide a realistic descrip- 
tion of the region of the solar atmosphere where the Fe 6173 A line (and the Ni 
6768 A line observed by MDI) are formed. The numerical simulation of the Sun's 
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Figure 1. Atlas profile of Fe 6173 A and six representative HMI filter transmission profiles. 



surface layers was computed with the radiation hydrodynamics code CO^BOLD 



(Freytag et a/., 2002 ( jWedemeyer et a/., 2004|). For a detailed documentation 



see http://www.astro.uu.se/~bf/co5bold_main.html, CO^BOLD solves the coupled 
non-linear equations of compressible hydrodynamics in an external gravity field 
together with the non-local radiative energy transport for a small fraction of the 
solar surface layers. A pre-tabulated equation of state, accounting for partial ion- 
ization and molecule formation, provides all necessary thermodynamic variables. 
The radiative transport is treated using realistic opacities adapted from a recent 



version of the MARCS stellar atmosphere package ('Plez et a/., 1992 Asplund et 



|a/., 1997 Gustafsson et a/., 2003|). T he simulation used in this study differs from 
the one used by Straus et al (2008) in two aspects: increased spatial resolution 
(both horizontal and vertical) and improved treatment of radiative transport. 
Compared to the simulation of [Straus et al. (2008), the present simulation has 



four times better horizontal resolution, three times better vertical resolution, 
and non-grey radiative transport treated in 12 frequency bands. To limit the 
amount of required core memory and time necessary to run the simulation, the 
horizontal box size was cut in half and the depth was also reduced. The simulated 
domain is defined by a fixed 3D Cartesian grid with 400 x 400 x 300 cells, each 
of size 14 km x 14 km x 7.5 km. The simulation thus covers a horizontal area 
of 5.6 X 5.6 Mm. Periodic lateral boundaries are used. The transmitting upper 
boundary is located at a height of approximately 900 km above the visible solar 
surface, while the open lower boundary lies 1.4 Mm below the surface. The part 
of the simulation that we use has completely relaxed and covers approximately 
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Table 1. Atomic parameters used for the line synthesis. 



element parameters"^ 



energy level parameters 



transition parameters 



log A 
W 

El. [eV] 
XI [eV] 
91. 
9\J 

Ao [A] 

log gf" 

Fr [rads-1] 

Fs [rads~^ cm^ 



Fe 6173 

7.52 
55.837 

2.2227 
7.870 

3 
1 

6173.3352 

-3.067 
2.042 X 10^ 
3.226 X 10-s 



Ni 6768 

6.30 
58.69 

1.826 
7.635 

1 

3 

6767.770 
-2.391 
1.0304 X 10^ 
8.276 X 10-^ 



" Partition functions are interpolated following 
ignored; ^ adapted values; ^ defined by Fvdw - 
the total hydrogen number density. 



Allen 



(19761; ^ Stark broadening is 



ITCt/WOO)^-^ Nnt, where A^Ht is 



two hours of solar time with snapshots saved every ten seconds. The code was 
executed on a 32-core node of an IBM-SP6 supercomputer and approximately 
three hours were necessary to advance ten seconds of solar time. 

A line-synthesis code, based on the assumption of local thermodynamic equi- 
librium (LTE) and written in IDL, was then fed with the physical parameters of 
a simulation snapshot to produce synthetic, two-dimensional spectra of several 
photospheric lines, including those of Fe I 6173 A and Ni I 6768 A. The atomic 
line parameters used are listed in Table 1. The oscillator strength and van der 
Waals coefficient were adapted in order to match the spatial mean spectrum 
to a solar atlas profile. Synthetic spectra were then calculated for all snapshots 
with this adapted parameter set. Using the recently published value of the iron 
abundance by Caffau et al (2011), we had to reduce the oscillator strength 
by approximately 35% with respect to the value given in the VALD database 
(Kupka et a/., 1999). Our description of the van der Waals broadening uses 



the modified exponent of 0.4 of the temperature dependence (Gomez et a/.. 



1987). The value of in Table 1 is approximately a factor two larger than 



one would get by translating the parameter given in Barklem et al (2000) to 



our reference temperature. While such modifications are not uncommon in line- 
synthesis calculations, it is interesting to note that they are still necessary in a 
3D simulation with such high spatial resolution. To understand the exact reason 
for this is not the goal of this work and difficult to resolve with only a single 
profile fit, but both direction and amount of the corrections seem compatible 
with the limits of the treatment in LTE and the uncertainties of the parameters. 

The main goal of this work is, instead, to determine the formation height of the 
HMI velocity signal in a realistic solar atmosphere, which is not a plane-parallel 
atmosphere but a highly corrugated surface. The line contribution functions, 
which were stored as a by-product of the line- synthesis calculations, give a quan- 
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titative idea of the size of this effect. Figure [2] shows the spatio-temporal averages 
of the intensity contribution functions of the hne cores (at the wavelength of 
minimum intensity) of Fe I 6173 A and Ni I 6768 A, and Figure [s] shows x- 
t diagrams of the height variations of the center-of-gravity of the contribution 
functions of these two lines. Figure [2] suggests that on average the cores of 
Fe 6173 and Ni 6768 are formed around 250 km and 300 km, respectively, in 
agreement with Bruls et al. |l990j and Norton et al (2006). However, as can 
be seen in Figure |3j there are large spatio-temporal variations of more than 
±150 km {i.e. several pressure scale heights). The formation height of the core 
of Fe 6173 varies between approximately 100 and 400 km, and that of Ni 6717 
between approximately 150 and 450 km. It is important to keep this in mind 
when talking about a "formation height". 
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3. Formation Height of the HMI Doppler Signal 

In the fohowing analysis we use only a ID (slit) cut through the 3D simula- 
tions. We multiplied the line profiles F(A,x,t) of the Fe I 6173 A line calculated 
in the model atmosphere with a representative set of HMI filter transmission 
profiles Ri{X^l < i < 6) and summed those products to construct filtergrams 
Ii{x,t;l < i < 6) = ^^F{X^x^t)Ri{X). In addition we determined the first 
moment (center of gravity) of the line profiles Vcg{x^t) = J^X F(A,x,t)dA / 
/^F(A,x,t)dA as well as the Doppler shift of the line core (Kore(^5^)) by de- 
termining the position of a fourth-order polynomial fitted to a range of ±32 mA 
around the minimum position. From the filtergrams Ii{x^t;l < i < 6) we 
determined Doppler velocities Vuuiix^t) using a simplified version of the HMI 
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Figure 3. x-t diagrams of the height variations (center-of- gravity) of the hne core intensity 
contribution functions of Fe I 6173 A (left) and Ni I 6768 A (right). The grayscale varies 
between 100 and 400 km for Fe 6173, and between 150 and 450 km for Ni 6768. 

observables-calculation pipeline, which comprises two steps: First, the phase of 
the first Fourier coefficient of the Fe I hne calculated from the six HMI intensities 
is computed, which yields an estimate of the Doppler velocity ("uncorrected 
velocity"): v = k x phase^ where /c is a constant. In a second step, "corrected" 
velocities are determined using a look-up table, which is based on the HMI 
filter-transmission profiles and an average Fe I 6173 A line profile (for details see 
[Couvidat a/., 2011) . 

Figure|4| shows an x-t diagram of the HMI Doppler velocities (left) and for 
comparison the actual vertical velocities in the simulation at a height of z = 
105 km (right). The corresponding x-t diagram of the center of gravity Doppler 
shifts Vcg{x^t) (not shown here) is practically indistinguishable from that of the 
HMI Doppler velocities. This good correspondence is also reflected by the very 
high correlation (Pearson linear correlation coefficient = 0.9991) between Vuui 
and Vcg. While the HMI Doppler shifts capture the spatio-temporal variations 
of the actual velocities very well, the contrast is significantly (24%) reduced: 
RMS(1/hmi) = 1195 ms-i vs RMS(l/io5km) = 1565ms-^ The RMS of Vcg is 
1254 ms~^, i.e. similarly reduced. 

Given the high correlation between Vhmi and Vcg, we believe that this re- 
duction is mainly an atmospheric modulation transfer function effect, which 
causes all perturbations that are not coherent over the height interval where 
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Figure 4. x-t diagrams of VuMli^, t) (left) and the actual velocities in the simulation at z = 
105km (right). Both have the same scaling (-4000ms-i < V < SOOOms-^). 



the signal is formed {e.g. high-frequency acoustic waves) to be smeared out and 
their amphtude reduced in the recovered Doppler signal. Another (although we 
believe minor) cause could be intrinsic calibration errors introduced through 
the correction with the look-up table. There are two issues: i) The RMS could 
be reduced because the slope of the look-up table is not correct and therefore, 
for a given phase of the first Fourier coefficient, the table returns too small a 
velocity. The slope could be incorrect because the average Fe I profile used to 
build the table is systematically different {i.e. broader and shallower) than the 
actual line profiles at a given pixel, ii) The look-up table is built using an average 
line profile and under the assumption that the only impact of a perturbation in 
the atmosphere is to shift the line in wavelength, without changing its shape. 
This, of course, is not the case, as the thermodynamic perturbations in the 
atmosphere do not just cause a Doppler shift, but also affect the width, depth, 
and asymmetry of the individual line profiles. 

To determine the formation height of the HMI Doppler signal, we calculated 
normalized cross-correlations CC{z) between ^hmi and the actual velocities 
V{z). The results are shown in Figurejsj together with the corresponding cross- 
correlation functions between Vcg and V{z)^ and Fcore and V{z). The cross- 
correlation function between Fhmi and V{z) is very similar to that between 
Vcg and V{z), conffiming the excellent correspondence between Fhmi and Vcg. 
This is perhaps not surprising, as the six HMI filters sample the whole line 
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profile, including the upper wings reaching into the continuum. Fhmi peaks near 
100 km, Vcg slightly below, in good agreement with t he velocity-response function 
calculated for the center of gravity of Fe 6173 A by Bello Gonzalez et al (2009). 
Therefore, one can conclude that the HMI Doppler signal is a very good measure 
for the center of gravity of the Fe 6173 A line, and is formed around 100km, i.e. 
rather low in the solar atmosphere. 

The cross-correlation function between Vcore and V{z) peaks near 230 km, 
slightly below the values given by Norton et a/. | (2006] ), but in good agreement 
with the NLTE value given by Bruls et al. (1990). The height of the maximum 
of the cross-correlation function between Kore and V{z) is sensitive to the wave- 
length range used for the polynomial fit. Extending that range from ±32 m A to 
±56 mA reduces the height of the maximum of the cross-correlation function to 
below 200 km. The reason for this is that by increasing the wavelength range the 
polynomials fitted extend into the line wings, which are formed deeper in the 
atmosphere. The highest layers in the atmosphere are probed only by a rather 
narrow wavelength range around the line core. 



4. Effects of Calibration Errors 

The HMI filter transmission profiles are believed to be known with an accuracy of 
the order of 1% (J. Schou, private communication). To obtain some insight into 
the effects of possible calibration errors on the measured velocities, we multiplied 
the blue (red) transmission profiles with a factor 1.01 (0.99), calculated new HMI 
intensities and from those determined new velocities using the same (unchanged) 
look-up tables as for the original velocities. That way we simulate unknown 
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changes of the transmission profiles that the pipehne process cannot account 
for. The same is repeated for assumed transmission profile errors of 0.3% and 
0.1%. The results are shown in Figure [g] 
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Figure 6. Differences between HMI Doppler shifts determined using the original filter trans- 
mission profiles and profiles that were modified by 1% (upper left), 0.3% (upper right), and 
0.1% (lower row). Only every fifth point is shown for clarity. 

The errors are not symmetric with respect to zero velocity and show a small 
velocity dependence. This is mainly due to the asymmetric error introduced, 
although the intrinsic asymmetry of the lines and the fact that the filters are 
not symmetric around zero wavelength might also play a role. For the case 
of a 1% (0.3%) [0.1%] change of the filter transmission profiles, the average 
offset introduced is about 170 ms~^ (60ms~^) [17ms~^], the RMS of the errors 
about 48ms~^ (17ms~^) [5ms~^]. We should emphasize, though, that these 
are worst case figures, since the actual HMI reduction pipeline comprises a third 
correction step, which makes use of the orbital velocity OBS_VR of SDO, which 
is known to great accuracy. The difference, as a function of OBS_VR, between 
this orbital velocity and the Doppler velocity averaged over the solar disk is 
fitted daily by a third-order polynomial. This polynomial is then used to correct 
the Dopplergrams and line-of-sight magnetograms. 

There are many possible sources of errors in the filter-transmission profiles. 
With above analysis only a single one is tested: a possible small error of the 
centering of the non-tunable part of the transmission profile of HMI, caused, 
e.g., by a de-centering of the blocking- filter transmission profile or an error in 
the phase of one of the Lyot elements. While a de-centering of the non-tunable 
part of the transmission profile is believed to be one of the largest possible error 
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Table 2. Random errors (in %) applied to the six HMI filter-response functions and 
resulting Doppler velocity offset and RMS errors (in ms"-*^). 





1 


2 


3 


4 


5 
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Offset 


RMS 
















[ms-l] 


[ms-l] 


Case 1 


0.76 


-0.59 


-0.37 


-0.025 


0.54 


0.96 
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Case 2 


-0.88 


-0.96 
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Figure 7. Differences between HMI Doppler shifts determined using the original filter trans- 
mission profiles and profiles that were modified by applying random errors of ±1% to the six 
HMI filters (c/. Table [I] for details). 



sources, there are others that appear more hke random errors. To gain some 
further insight, we have therefore studied the effects of randomly distributed 
errors by multiplying the six HMI filters with factors (1 + 6^), with Ci {1 < i < 6) 
being random numbers between -0.01 and +0.01 {cf. Table [4|. The results are 
shown Figure [7| The offset for the three (random) cases varies between -65 
and +35ms~^ and the RMS error between 10 and 30ms~^. The slopes of the 
velocity errors are very different for the three cases, and even change sign. This 
demonstrates that the asymmetry of the errors introduced for the tests leading 
to Figure [6] is more important than the intrinsic asymmetry of the Fe line. 
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Figure 8. Atlas profile of Ni 6768 A and five representative MDI filter transmission profiles, 
four of which are used for the velocity determination and one (in red) to measure the continuum 
intensity. 



5. Formation Height of the MDI Doppler Signal 

The HMI instrument on SDO replaced the MDI instrument on SOHO, which 
has recorded solar oscillations for over an entire solar cycle. The HMI instrument 
observes in a different spectral line than MDI and uses more and narrower filter 
settings. Therefore, it is reasonable to expect that the HMI and MDI Doppler 
signals form at different heights in the solar atmosphere. Since the HMI and 
MDI series might be combined for solar-cycle related studies, it is of interest 
to know how large this difference is. 'Yurchyshyn and Wang ( 2QQ1| adopted a 
formation height of 200 km for the MDI Doppler signal based on a study by Jones] 
(|1989|), [Straus et aL\^00^ assumed 100 km, [Wachter| (|2008|) in his study on the 



instrumental response function of filtergraph instruments obtained 180 km for 
the MDI signal for the quiet-Sun reference model of Maltby et al. (1986), and 
others simply use the Ni 6768 A core formation height of 300 km. To clarify this 
issue and get a better idea of the formation height of the MDI Doppler signal, 
we repeated the procedure as described in Section [s] for the Ni I 6768 A line 
and the theoretical MDI filter-response functions, which are shown in Figure [8] 
(R. Wachter, private communication). 

Figure [9] shows the cross-correlation between the Doppler shifts Fmdi "mea- 
sured" by MDI and the actual vertical velocities in the model atmosphere [V^(^)], 
together with the corresponding cross-correlation functions between V^g of Ni 
6768 A and V{z)^ and Kore of Ni 6768 A and V{z). Comparing this diagram to 
Figure |5] one can note several differences: i) the cross-correlation of the MDI 
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Doppler signal (Fmdi) is markedly different from that of the center of gravity of 
Ni 6768 A (Kg); this is due to the fact that MDI uses only four filter positions and 
therefore does not sample the full profile of Ni 6768 A; ii) the cross-correlation 
of Vmbi has its maximum at about 125km, i.e. about 25km higher than that 
of "i^Mi; this is slightly higher than the value used by Straus et al. (2008), but 
55 and 75 km lower than the values in Wachter (2008) and Yurchyshyn and 



Wang (2001); and in) the cross-correlation of Kore peaks at around 300 km, 



significantly higher than the corresponding function of Fe 6173, and in excellent 
agreement with Norton et al ( |2006, ), who quoted 288 km for the core formation 
height of the Ni line. 



6. Effects of Spatial Resolution 

With a horizontal grid size of 14 km, the 3D simulations that we used have 
a much higher spatial resolution than the actual HMI (or MDI) observations. 
Small-scale features that are invisible in HMI observations but prominent in 
the simulations at certain heights might skew the results. To investigate the 
effects that limited spatial resolution may have on the formation height of the 
HMI and MDI Doppler signals, we repeated the analysis presented in Sections 
|3] and |5] after convolving the Fe and Ni spectra in the spatial dimension with a 
Gaussian of 1.5^' full width at half maximum (FHWM). This figure is based on a 
simple empirical model of the HMI point spread function (PSF) that was fitted 
to HMI calibration data ( [Wachter et a/., 2011| ). The model includes a narrow 
core (described by a Gaussian) and an extended tail (accounting for scattered 



light). Wachter et al (2011) give a 1/e width of the central Gaussian core of 1.8 



pixel, equivalent to a FWHM of approximately 1.5^' 
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Figure 10 shows the resulting x-t diagrams for the HMI (left) and MDI (right) 



Doppler signals, Figure 11 the corresponding cross-correlation functions with the 
spatially convolved vertical velocities in the atmosphere. 



100 




2000 



1000 



hi 000 



L-2000 



x/ Mm 



2 3 4 
x/ Mm 



Figure 10. x-t diagrams of ViiMi(x,t) (left) and V-MDiix^t) (right) after convolving the 
Fe and Ni spectra with a Gaussian of 1.5" FWHM. The scaling is different from that in 
Figure |4] To account for the significantly reduced RMS, the range has been reduced to 
(-2000ms-i < V < 2000ms-i). 

The two x-t diagrams of the HMI and MDI Doppler signals in Figure [lO] 
look very similar, as expected given the small difference in formation height of 
the two signals. Comparing them to the full-resolution x-t diagram in Figure [4j 
one can notice striking differences. Most of the small-scale convective motions 
are lost in the signal when reduced to actual HMI resolution and the spatio- 
temporal pattern becomes dominated by the larger-scale five- minute p modes. 
As a result of this different weighting of convective motions and p modes, the 
apparent formation height increases by about 50 km to 150 km for the HMI signal 
and to 175 km for the MDI signal (see the positions of the maxima of the solid 
lines in Figure 11 ). 



We should emphasize that the model PSF of Wachter et al. (2011) is based 



on pre-launch data, not on measurements in space. Several factors typical for 
ground measurements (air flows, temperature changes, etc.) may have impacted 
the image quality. A PSF of 1.5'^ is therefore a worst case, and it is likely that the 
actual HMI PSF is better. To quantify the range of the effects of limited spatial 
resolution, we repeated the above procedure assuming the PSF of a diffraction- 
limited 14-cm telescope operating at 6173 A, which can be approximated by a 
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Figure 11. Left: Cross-correlations between Vumi, V'cg, V'core and the actual velocities V(z) 
as a function of height in the model atmosphere after convolving the Fe 6173 spectra and 
vertical velocities with a Gaussian of 1.5" FWHM. Right: Same for the MDI Ni 6768 A line. 

Gaussian of 0.9'' FWHM. In that case the apparent formation height of both 
the HMI and MDI Doppler signal is found to be about 40 km higher than in 
the fuh-resolution case. With the actual HMI PSF falling somewhere in between 
the 0.9'' and 1.5" case, one can therefore assume a height effect due to spatial 
smearing of 40 to 50 km. 

7. Conclusions 

We have studied the formation height of the HMI and MDI Doppler signals using 
3D radiation-hydrodynamic simulations. From the synthesized Fe 6173 and Ni 
6768 A lines HMI and MDI Doppler signals were calculated using representative 
filter-response functions and (simplified) pipeline processes. The resulting HMI 
Doppler signal is nearly identical to the center-of-gravity shift of the Fe 6173 A 
line (linear correlation coefficient = 0.9991). It is formed at a height of about 
100 km, about 25 km lower than the MDI Doppler signal and significantly lower 
than some previous estimates. The formation height shows a small dependence 
on the spatial resolution and increases by 40 - 50 km to about 140 - 150 for the 
HMI signal and to about 165 - 175 km for the MDI signal. An error of 0.3% 
in the HMI filter transmission profiles (probably more realistic than a 1% error, 
thanks to the a-posteriori polynomial correction applied to some observables) 
leads to an error of less than 20ms~^ (RMS) in the measured HMI Doppler 
velocities. The current analysis is applicable only to disk center (/i = 1) observa- 
tions. The inclusion of center-to-limb effects would significantly complicate the 
foregoing analysis. A more interesting study would be to expand this analysis 
from HD (hydrodynamics) to MHD (magnetohydrodynamics) and investigate 
the magnetic-field diagnostics potential of HMI. 

Acknowledgements The simulations were carried out at CINECA (Bologna/Italy) with 
CPU time assigned under INAF/CINECA agreement 2008/2010. SDO is part of NASA's Living 
With a Star (LWS) program. HMI was designed and assembled at Stanford University and 
Lockheed Martin Solar and Astrophysics Laboratory. We are grateful to the SDO/HMI team 
for making this instrument a reality. S.C. was supported by NASA Grant NAS5-02139 (HMI). 
Th.S. acknowledges financial support by ASI. SOHO is a mission of international cooperation 



accepted for publication in Solar Physics: 27 April 2011; p. 14 



On the Formation Height of the SDO/HMI Doppler Signal 



between ESA and NASA. We thank G. Severino for helpful discussions and an anonymous 
referee for constructive comments, which helped to improve the paper. 



References 

Allen, C. W. 1976, Astrophysical Quantities, University of London, The Athlone Press, 3rd 
edition. 

Asplund, M., Gustafsson, B., Kiselman, D., Eriksson, K. 1997, Astron. Astrophys., 318, 521. 
Barklem, P.S., Piskunov, N., O'Mara, B.J.: 2000, Astron. Astrophys. Suppl. 142, 467. 
Baur, T.G., House, L.L., Hull, H.K.: 1980, Solar Phys. 65, 111. 
Bello Gonzalez, N., Okunev, O., Kneer, P.: 1996, Astron. Astrophys. 490, L23. 
Bello Gonzalez, N., Yelles Chaouche, L., Okunev, O., Kneer, P.: 2009, Astron. Astrophys. 
494, 1091. 

Bruls, J.H.M.J., Lites, B.W., Murphy, G.A.: 1991, In: November, L.J., (ed.). Solar Polarimetry, 

National Solar Observatory, Sunspot, NM, 444. 
Caffau, E., Ludwig, H.-G., Steffen, M., Preytag, B., Bonifacio, P.: 2011, Solar Phys. 268, 255. 
Couvidat, S., Schou, J., Shine, R.A., Bush, R.I., Miles, J.W., Scherrer, P.H., Rairden, R.L.: 

2011, Solar Phys., DOI: 10.1007/sll207-011-9723-8. 
Domingo, V., Pleck, B., Poland, A.I.: 1996, Solar Phys. 162, 1. 
Preytag, B., Steffen, M., Dorch, B. 2002, Astrono. Nach. 323, 213. 

Gomez, M.T., Severino, G., Marmolino, C., Roberti, G.: 1987, Astron. Astrophys. 188, 169. 

Gustafsson, B., Edvardsson, B., Eriksson, K., Mizuno-Wiedner, M., J0rgensen, U.G., Plez, 
B.: 2003, IniHubeny, I., Mihalas, D., Werner, K. (eds).. Astronomical Stellar Atmosphere 
Modeling, CS-288, Astron. Soc. of the Pacific, San Francisco, 331. 

Jones, H.P.: 1989, Solar Phys. 120, 211. 

Kupka, P., Piskunov, N., Ryabchikova, T.A., Stempels, H.C., Weiss, W.W.: 1999, Astron. 

Astrophys. Suppl. 138, 119. 
Maltby, P., Avrett, E.H., Carlsson, M., Kjeldseth-Moe, O., Kurucz, R.L., Loeser, R.: 1986 

Astrophys. J. 306, 284. 

Muller, R., Dollfus, A., Montague, M., Moity, J., Vigneau, J.: 2000, Astron. Astrophys. 359, 
373. 

Norton, A. A., Pietarila Graham, J., Ulrich, R.K., Schou, J., Tomczyk, S., Liu, Y., Lites, B.W., 
L'opez Ariste, A., Busch, R.L, Socas-Navarro, H., Scherrer, P.H.: 2006, Solar Phys. 239, 
69. 

Plez, B., Brett, J.M., Nordlund, A.: 1992, Astron. Astrophys., 256, 551. 

Scherrer, P.H., Bogart, R.S., Busch, R.L, Hoeksema, J.T., Kosovichev, A.G., Schou, J., Rosen- 
berg, W., Springer, L., Tarbell, T.D., Title, A., Wilfson, C.J., Zayer, L, and the MDI 
engineering team: 1995, Solar Phys. 162, 129. 

Schou, J., Scherrer, P.H., Bush, R.L, Wachter, R., Couvidat, S., Rabello-Soares, M.C., et al: 
2011, Solar Phys., in preparation. 

Straus, T., Pleck, B., Jefferies, S.M., Cauzzi, G., Mcintosh, S.W., Reardon, K., Severino, G., 
Steffen, M.: 2008, Astrophys. J. 681, L125. 

Wachter, R.: 2008, Solar Phys. 251, 491. 

Wachter, R., Schou, J., Rabello-Soares, M.C., Miles, J.W., Duvall, T.L., Bush, R.L: 2011, Solar 

Phys., DOI: 10.1007/sll207-011-9709-6. 
Wedemeyer, S., Preytag, B., Steffen, M., Ludwig, H.-G., Holweger, H. 2004, Astron. 

Astrophys., 414, 112. 
Yurchyshyn, V.B., Wang, H.: 2001, Solar Phys. 202, 309. 



accepted for publication in Solar Physics: 27 April 2011; p. 15 



accepted for publication in Solar Physics: 27 April 2011; 



